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Abstract 

We summarize the results of a recent global analysis of proton and 
deuteron F<i structure function world data performed over a large range 
of kinematics, including recent measurements done at JLab with the CLAS 
detector. From these data the lowest moments (n < 10) of the unpolarized 
structure functions are determined with good statistics and systematics. 
The Q 2 evolution of the extracted moments is analyzed in terms of an OPE 
based twist expansion, taking into account soft-gluon effects at large x. A 
clean separation among the Leading and Higher- Twist terms is achieved. 
By combining proton and deuteron measurements the lowest moments of 
the neutron F2 structure function are determined and its leading twist term 
is extracted. Particular attention is paid to nuclear effects in the deuteron, 
which become increasingly important for the higher moments. Our results 
for the non-singlet, isovector (p — n) combination of the leading twist mo- 
ments are used to test recent lattice simulations. We also determine the 
lowest few moments of the higher twist contributions, and find these to be 
approximately isospin independent, suggesting the possible dominance of ud 
correlations over uu and dd in the nucleon. 



a To appear in the Proceedings of the TV International Conference on Quark and Nuclear 
Physics (QNP06), Madrid (Spain), June 5-10, 2006. 



In the past few years, thanks to the performance of the CLAS detector of the 
Hall B at Jefferson Lab, the unpolarized structure functions F2(x, Q 2 ) of both pro- 
ton and deuteron have been measured precisely in a wide continuous interval of 
values of the Bjorken variable x and of the squared four-momentum transfer Q 2 
[HI2|- In particular, the F 2 structure functions have been extracted over the whole 
resonance region (W < 2.5 GeV) below Q 2 ~ 4.5 (GeV/c) 2 . These measurements, 
together with existing world data, have allowed the determination of the moments 
M n (Q 2 ) of the structure function F 2 up to n = 10, drastically reducing the un- 
certainties related to data interpolation and extrapolation as well as providing 
the most accurate evaluation of the Q 2 -dependence of the moments starting from 
Q 2 ~ 0.1 (0.5) (GeV/C) 2 for proton (deuteron) up to Q 2 ~ 100 (GeV/c) 2 . 

A very powerful tool for analyzing in QCD the Q 2 -dependence of the moments 
M n (Q 2 ) is represented by the Operator Product Expansion (OPE). Through OPE 
the moments are given as a series expansion in terms of matrix elements of local 
operators, namely 



where /i is the renormalization scale, O nT (/i) is the reduced matrix element of 
local operators with definite spin n and twist r, related to the non-perturbative 
structure of the target, and E nr (fi,Q 2 ) is a short-distance coefficient calculable in 
pQCD thanks to the asymptotic freedom. 

The first term of the series corresponds to the Leading Twist (LT) and 
is determined by single-parton distributions in the target. Subsequent terms are 
the so-called Higher Twists (HT's), which depend on the interactions among par- 
tons (i.e., on multi-parton correlations). Their determination is considerably more 
challenging both experimentally and theoretically. 

In Refs. P and [2] both the LT and the HT's (up to twist-6) have been extracted 
from the proton and deuteron moments, respectively. Two important features of 
our twist analysis should be mentioned: 1) the analyzed moments are Nachtmann 
moments, and not Cornwall-Norton ones; 2) the perturbative evolution of the LT 
is calculated using soft-gluon resummation techniques. 

We have adopted the Nachtmann definition of the moments in order to get rid of 
the target-mass corrections, which are power corrections of kinematical origin and 
affect the Cornwall- Norton (CN) definition of the moments. The use of Nachtmann 
moments allows to extract in a clean way dynamical HT's, which are the only 
ones related to multi-parton correlations. The target-mass corrections can be 
subtracted in principle from the CN moments, but this procedure requires a precise 
knowledge of the LT moments at large n. Moreover, the contribution of target- 
mass corrections is sizable just in the region of Q 2 ~ few (GeV/c) 2 , where HT's 



M n (Q 2 ) = Y / E nT (fx,Q 2 ) OnM 




(1) 



r=2 



2 



are expected to show up. This point is clearly illustrated in Fig. ^ where the 
Nachtmann moment m\ D \q 2 ) is compared with the corresponding CN one as 
determined in Ref. It is also worthwhile to note the much larger number of 
very accurate data points obtained thanks to the CLAS detector. 
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Figure 1: Moment Mg (Q 2 ) of the deuteron versus Q 2 . Open dots are the results 
of Ref. |1J/ using the Nachtmann definition of the moments, while full squares 
correspond to the CN moment as determined in Ref. J^j. The reported systematic 
errors have been precisely estimated only in Ref. Jjfl. 

As for the perturbative evolution of the LT, the presence of large logarithms in 
the coefficient function at large n requires to go beyond any fixed-order approxima- 
tion. As shown in Ref. j^j for a reliable extraction of both LT and HT's it is crucial 
to apply the resummation of soft gluons. Moreover, for internal consistency the 
value of a s (M§) has been extracted in Ref. |B] by analyzing the large Q 2 behavior 
of the proton moments of Ref. £Q at the next-to-leading logarithmic accuracy used 
for the soft-gluon resummation. 

The extracted LT components of proton and deuteron moments can be com- 
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bined to form the LT moments of the neutron F 2 structure function. In the Impulse 
Approximation the deuteron structure function can be written as the sum of two 
terms: = p( D ' conv ^ _|_ ^W/^ ^gj-g the fi rs t term is given by a convolution 

of the nucleon structure function and the nucleon momentum distribution func- 
tion in the deuteron, f^ D \ while the second term represents nuclear off-shell plus 
relativistic corrections 13 . Thus, the neutron moments can be obtained from: 

1 _ A(o//) 

Mi n \Q 2 ) = 2Mf)(Q 2 ) " - M^(Q 2 ) , (2) 

M ' 

where is the moment of the function and A^°-^ represents the nuclear 
off-shell correction. Thus the quantity /[l — is the global nuclear 

correction factor, i.e. the EMC effect of the deuteron in moment space. The 
nuclear corrections and their uncertainties have been carefully estimated in Ref. [Sj. 

The extracted neutron LT moments can be combined with the proton LT ones 
to form the non-singlet, isovector (p — n) LT moments of the nucleon F 2 structure 
function. Our findings are compared against lattice QCD simulations in Fig. |21 
While a linear extrapolation of the lattice results to the physical pion mass overes- 
timates our data significantly, the results obtained using the extrapolation method 
of Ref. [H] are in much better agreement. Our results for higher moments exhibit 
a remarkable precision, and it would be valuable to compare them with corre- 
sponding lattice moments, because the effects of chiral loops are expected to be 
suppressed at large n. 

Another interesting issue is the behavior of the neutron to proton structure 
function ratio R(x) = F^ix) / F 2 {x) in the limit x — > 1. This ratio can be easily 
related to the corresponding d/u ratio of parton distributions. The "standard" 
value R(x = 1) = 1/4 is used in most of the parton distribution fits and corresponds 
to a vanishing d/u ratio, while the value R(x = 1) = 3/7 is expected from pQCD 
arguments and corresponds to a limiting value of 1/5 for the d/u ratio. As shown 
in Ref. [Sj the ratio of neutron to proton moments, M™/M%, is largely independent 
of Q 2 and is very sensitive at large n to the large-x behavior of R(x). Our results 
are shown in Fig. El The trend is clearly towards the "standard" value 1/4 as 
n increases, although the precision of the data does not exclude higher limiting 
values. From Fig. El one can also see the impact of the nuclear corrections on the 
deuteron moments at large n: for n = 12, which corresponds to an average value 
of x around 0.75, the nuclear correction corresponds to a reduction factor of ~ 2, 
though within large uncertainties. 

Once the LT contribution to moments is determined, one can study the isospin 
dependence of the HT contribution. Assuming that final state interactions and me- 

b For the moment with n = 2 shadowing and meson-exchange-current should be also considered 
(see Ref. 0). 
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Figure 2: Isovector (p — n) moments of the nucleon F 2 structure function compared 
with lattice QCD simulations at Q 2 = 4 (GeV/c) 2 . Full dots are the data from 
Ref- open diamonds and triangles correspond to the lattice results from Refs. ^ 
and Jjfl, obtained assuming a linear extrapolation in the quark masses, and open 
squares are the results of Ref. where chiral loop effects are considered. 



son exchange currents have negligible impact on the HT's above Q 2 ~ 1 (GeV/c) 2 , 
the same nuclear corrections used for the LT can be applied to the HT's. In this 
way the total HT contributions to proton moments and to the corrected deuteron 
moments can be compared, as done in Fig. |3J The comparison indicates clearly 
that the total HT contribution is almost independent of isospin, suggesting the 
possible dominance of ud correlations over uu and dd in the nucleon More- 
over, the isovector combination (p — n) of the F2 structure functions should be 
almost free from HT contributions within the present uncertainties. 

In conclusion we have analyzed experimental data on proton and deuteron 
F 2 structure functions in order to determine their moments. Using the OPE we 
have extracted both leading and higher twist contributions. By combining proton 
and deuteron moments we have determined the LT moments of the neutron F2 
structure function, paying particular attention to the issue of nuclear effects in the 
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deuteron, which are increasingly important at large n. 

The main results of our analysis can be summarized as follows: 

• the ratio of neutron to proton moments is consistent with F^jF^ — > 1/4 as 
x — ► 1, although the higher value of 3/7, suggested by pQCD arguments, 
cannot be excluded; 

• the non-singlet moments are in good agreement with the lattice data |H] 
only when the latter are extrapolated to physical quark masses after having 
taken into account chiral loops P; 

• the total contribution of HT's is found to be isospin independent, suggesting 
the possible dominance of ud correlations over uu and dd in the nucleon. This 
implies that in the isovector combination (p — n) of F 2 structure functions 
the HT's are almost consistent with zero within the present uncertainties. 
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Figure 3: Ratio of neutron to proton moments versus n: full triangles and open 
dots corresponds to the results of Ref. J3J/ obtained with and without the nuclear 
corrections, respectively. The solid (dashed) line indicates the scenario where 
d/u -> (1/5) forx -»• 1. 
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Figure 4: Total higher twist contribution to the proton (filled triangles) and 
deuteron (open squares) moments evaluated at Q 2 = 2 (GeV/c) 2 . 
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